CONSTRUCTION OF POLYNOMIAL LATTICE RULES OVER F 2 
WITH SMALL MEAN SQUARE WEIGHTED C 2 DISCREPANCY* 

TAKASHI GODAt 

Abstract. The C2 discrepancy is one of the well-known quantitative measures for the irregularity 
of distribution of point sets in the unit cube. The weights were introduced by Sloan and Wozniakowski 
to take into account the relative importance of the discrepancy of lower dimensional projections. As 
known under the name of quasi-Monte Carlo methods, point sets with small weighted C2 discrepancy 
arc of use in numerical integration. In this study, we study the component- by-component construction 
of polynomial lattice rules over the finite field F2 whose scrambled point sets have small mean square 
weighted £2 discrepancy. We prove an upper bound on this discrepancy which converges at almost 
the best possible rate of N~ 2+s for all <5 > 0, where N denotes the number of points. Numerical 
experiments confirm that the performance of our constructed polynomial lattice rules is comparable 
or even superior to that of the well-known Sobol' sequences. 
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1. Introduction. In this paper, we study the approximation of an s-dimensional 
integral over the unite cube [0, l) s 

/(/) = / f(x)dx, 

by averaging with equal weights the function values evaluated at N points 

1 £z? 

Q(f) = 

n=0 

Monte Carlo (MC) and quasi-Monte Carlo (QMC) methods choose the point set 
Pn.s = {^o, ■ ■ ■ , xn-i} randomly and detcrministically, respectively. The aim of 
QMC methods is to distribute the quadrature points as uniformly as possible so as 
to yield the smaller integration error. This idea is supported by the general form of 
various integration error bounds 

\I(f)-Q(.f)\<V(f)D(P N , s ), (1.1) 

where V(f) is the variation of the integrand / in some sense, which depends only 
on /, while D(Pn. s ) is the corresponding discrepancy of the point set Pn,s, which 
measures the irregularity of distribution of Pn.s and depends only on Pn.s- Thereby, 
as the smaller D(Pn. s ) is, the smaller integration error we can expect. The most well- 
known bound is given by the so-called Koksma-Hlawka inequality in which V(f) is the 
variation of / in the sense of Hardy and Krause and D(Pn, s ) is the star discrepancy 
of Pn.s, see for example [T71 120] . 

Randomization of the QMC point set is helpful to obtain statistical information 
on the integration error and sometimes even enables us to improve the rate of con- 
vergence for numerical integration. There have been several methods introduced for 
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randomization [5J [T5J [25] 13S] ■ Using the linearity of expectation and the mean 

square integration error is upper-bounded by 



where the expectation is taken with respect to all the possible randomized point sets 
of Pn,s- Hence, the mean square discrepancy becomes a meaningful measure of the 
irregularity of distribution of -Pj\r,s in this setting. 

Among others for D(Pn, s ), the Li discrepancy is one of the popular measures 
of the irregularity of distribution of point sets. The relationship between the Li 
discrepancy and numerical integration has been often discussed in the literature, see 
for example [HI EH EH [38] [39] . Sloan and Wozniakowski [34] introduced a concept 
of the weighted Li discrepancy to take the relative importance of the discrepancy 
of lower dimensional projections into account. It provides a part of the reason why 
QMC methods are successful even for very large values of s, as often reported in the 
practical applications to financial problems 3, 22, 29 . This phenomenon is hard to 
explain by the classical integration error bounds. Hence, construction of point sets 
with small weighted Li discrepancy is of particular interest to QMC practitioners. 
Especially, in this paper, we focus on constructing randomized QMC point sets with 
small mean square weighted Li discrepancy. 

In order to give the definition of the weighted Li discrepancy, we introduce some 
notation first. For a point set Pn.s = {xq, ■ ■ ■ , &n-i} in the unit cube [0, the 
local discrepancy function is defined as 



where t = (ii, . . . ,t s ) is the vector from [0, l) s , [0,t) is the axis-parallel box of the 
form [0,ti) x • • • x [0, t s ), and An([0, i), Pjv.s) denotes the number of indices n with 
x n £ [0, t). Let I s = {1, . . . , s} and let 7„ be a non-negative real number for u C I s . 
We denote by |u| the cardinality of u and by t u the vector from [0, l)' u containing all 
the components of t £ [0, l) s whose indices are in u. Further, let dt u = Ylj eu dtj and 
let (t u , 1) denote the vector from [0, l) s with all the components whose indices are 
not in u replaced by one. Then, the weighted Li discrepancy of the point set Pn,s is 
given by 



We can recover the classical Li discrepancy by choosing jj a = 1 and 7„ = for 
u C I s - The following proposition generalizes the well-known formula for the classical 
Li discrepancy introduced by Warnock, see for example [H [19] . 

Proposition 1.1. For any point set Pn.s = {xo> ■ • • > fjv-i} in [0, l) s and any 

sequence 7 = (7„) u cj s of weights, we have 

^■2,n,-t(Pn,s) 



E [|/(/) - Q(/)| 2 ] < V 2 (f)E [P> 2 (Pn, s )] , 





1 2 

3~H ~ N 




+ 



1 

iV2 



n,n'— j&u 



N-1 



JJ(1 -max(x ni j,x n >j)) 



where x n _j is the j-th component of the point x 
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There are two prominent construction principles of the QMC point sets: lattice 
rules [101131] and digital (i, m, s)-nets [SI HO]- In this study, we are concerned with 
polynomial lattice rules which can be categorized into the latter, while its name 
comes from the analogy with lattice rules. Since first introduced by Niederreiter [21 j . 
polynomial lattice rules have been extensively investigated, sec for example [8, 16 . In 
the following, we give the definition of polynomial lattice rules for the case of base 2 
because we will only deal with that case. 

Let F2 := {0, 1} be the two element field and represent by ¥2((x~ 1 )) the field of 
formal Laurent series over F2. Every element of F2((x -1 )) has the form 



l—w 



where w is an arbitrary integer and all t; G F2. Further, we denote by ¥2[x] the set of 
all polynomials over F2. For a given integer m, we define the map v m from F2((a; -1 )) 
to the interval [0,1) by 

(oc \ rn 

l—w / i=max(l,Kj) 

We often identify a non-negative integer k whose dyadic expansion is given by k = 
kq + k±2 + ■ ■ - + ^02° with the polynomial over F%[x] as k(x) = k,$ + i%ix ■ ■ - + K a x a . For 
k = (ki ■ ■ ■ , k s ) G F 2 [:e] s and q = (qi ■ ■ ■ , q s ) G F2[x] s , we define the inner product as 

s 

k ■ q = kjqj G F 2 [x], 
i=i 

and we write q = (mod p) if p divides q in F 2 [a;]. Using these notations, the 
polynomial lattice point set is constructed as follows. 

Definition 1.2. Let m,s G N. Let p G F2[x] be an irreducible polynomial 
with deg(p) = m and let q = (gi, . . . , q 3 ) G F2[x] s . The polynomial lattice point set 
P2 m ,s{q,p) is the point set consisting of the 2 m points 

( ( n{x)q 1 {x) \ ( n(x)q s {x) \\ 

for n G ¥2[x] with deg(n) < m. 

In the following, the notation P2 m ,s(?,p) implicitly means that deg(p) = m and 
the number of components for a vector q is s. 

For randomization of the polynomial lattice point set, we consider to apply Owen's 
scrambling [26J [27l [28]. It proceeds as follows. For x = (xi,--- ,x s ) G [0, l) s , we 
denote the dyadic expansion by xj = x.j.xlr 1 + .Tj. 2 2 _2 + • • • . Let y = (yi, ■ ■ ■ , y 8 ) G 
[0, l) s be the scrambled point of x whose dyadic expansion is represented by yj = 
yj \2r x + ?/j j2 2~ 2 + • • • . Each coordinate yj is obtained by applying permutations to 
each digit of Xj. Here the permutation applied to xj^ depends on xjj for 1 < I < k— 1. 
In particular, y jtl = Kj(xj t i), y jt 2 = ^j,x jtl (^,2), Vj,3 = ^j.xj, 1,^,2 (^,3), and in general 

Vj,k = n j,x jt i,--- ,Xj ik -! (Xj,k)i 

where itj^j 1,— ,x ■ is a random permutation of {0,1}. We choose permutations 
with different indices mutually independent from each other where each permutation 
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is chosen with the same probability. Then, as shown in [261 Proposition 2], the 
scrambled point y is uniformly distributed in [0, l) s . 

Our aim here is to find a vector q with p fixed, which yields small mean square 
weighted £2 discrepancy Restricting each qj G F2[a;] such that deg (</_,-) < m and qj ^ 
0, the number of candidates for q is (2 m — l) s , which is quite large. The componcnt- 
by-componcnt (CBC) construction can significantly reduce the computational burden 
by searching over all the candidates of qj+i while leaving the existing components 
{qi, ■ ■ ■ ,Qj) unchanged. The CBC construction was first invented for lattice rules 
by Korobov [12] and re-discovered more recently by Sloan and Reztsov [35]. It also 
has been applied to polynomial lattice rules. Without requiring exhaustive search, 
the CBC construction usually finds a good vector q as discussed in many previous 
studies, see for example [5] [5] [TU [TS] . Hence, we employ the CBC construction to 
find a vector q which gives small mean square weighted £2 discrepancy. 

We end this section with a brief outline of this paper. In the next section, we 
introduce Walsh functions and their useful properties. They play a central role in 
the analysis of the mean square weighted £2 discrepancy. In Section |3j we study the 
mean square weighted £2 discrepancy of scrambled polynomial lattice rules. Next, 
in Section [4] we construct polynomial lattice rules which have small mean square 
weighted £2 discrepancy. We prove an upper bound on the root of this discrepancy 
which converges at a rate of N~ 1+s for all 6 > 0, where N = 2 m denotes the number 
of points. As Roth [31] and Proinov [30] proved that the lower bound on the classical 
£2 discrepancy of N points is given by 

(l0gJV)(«-l)/2 

£2,jv, T (-rN>) > c s — , (1.2) 

where c s is a constant dependent only on s, our upper bound is almost best possible 
in the sense that a rate of ./V -1 cannot be achievable. We further mention about the 
strong tractability of our construction algorithm. Finally, in Section [5] we show the 
performance of our constructed polynomial lattice rules and compare with that of the 
well-known Sobol' sequences. 

2. Walsh functions. Walsh functions were first introduced by Walsh [37] and 
have been extensively studied for example in [4j [10] . We refer to [8] Appendix A] for 
more information on Walsh functions. In the following, No := NU{0} denotes a set of 
non-negative integers. We first give the definition of dyadic Walsh functions for the 
one-dimensional case. 

Definition 2.1. Let k e No with dyadic expansion k = n + k\2 + ■ ■ ■ + n a 2 a . 
Then, the k-th dyadic Walsh function wah : [0, 1) — > { — 1, 1} is defined as 

walfc(x) = (_i)*i«o+-+*-+i«« J 

for x G [0, 1) with dyadic expansion x = xil^ 1 + iE22~ 2 + • • • (unique in the sense 
that infinitely many of the Xi must be zero). 

This definition can be generalized to higher-dimensional case. 

Definition 2.2. For s e N, let x = (xi, ■ ■ • ,x s ) e [0, l) s and k = (ki, ■ • ■ , k s ) 6 
Ng. We define wal fc : [0, l) s {-1, 1} by 



walfe(a;) = wal fcj (xj). 
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In the following, the operator © denotes the digitwisc addition modulo 2, that is, 
for x,y G [0, 1) with dyadic representations x — YliLi x i^~ % an d y — Vi^~' ''■> © 
is defined as 

oo 

x®y = ^z{T\ 

1=1 

where Z{ = Xi + yi (mod 2). We also define a digitwise addition for non- negative 
integers based on those dyadic representations. In case of vectors in [0, l) s or N§, 
the operator © is carried out componentwise. Further, we call x G [0, 1) a dyadic 
rational if it can be represented by a finite dyadic expansion. The proposition below 
summarizes some basic properties that Walsh functions hold. 
PROPOSITION 2.3. We have the following : 

1. For all k,l G N and all x, y G [0,1) with the restriction that if x,y are not 
dyadic rationals, then x © y is not allowed to be a dyadic rational, we have 

w&\ k {x)w&\i(x) = wal fce i(a;), wal fc (x)wal fc (?/) = wah(x © y). 

2. We have 

i f i 

waio(x)dx = 1 and / walk(x)dx = if k £ N. 



1 ifk = l, 
erwise. 



3. For all k,l € Ng ; we have 



4- For s G N, the system {walfc : k = (fci . . . , k s ) G Nq} is a complete orthonor- 
mal system in ^([O, l] s )- 
Furthermore, in order to introduce an important relation between Walsh functions 
and polynomial lattice rules as described below in Lemma 12.51 we add one more 
notation and introduce a concept of the so-called dual net of a polynomial lattice 
point set P2™,s(QiP)- Fot k G No with dyadic expansion k = k + k\2 + • • ■ , tr m (fc) 
gives a polynomial of degree at most m by truncating the associated polynomial 
k(x) G ¥ 2 [x] as 

tr m (fc) = k Q + kix + • • • + k m -ix m ~ . 

For a vector k = (ki ■ ■ ■ ,k s ) G Nq, we define tr m (fe) = ( ti m (ki), • • ■ , tr m (k s )). With 
this notation, we introduce the following definition of the dual net D* p . 

Definition 2.4. The dual net for a polynomial lattice point set P2^. s (q,p) is 
given by 

D* q p = {k E N s : tr m (fc)-q = (mod p)}. 

Then, the following lemma relates the dual net of a polynomial lattice point set to 
the numerical integration of Walsh functions. It follows immediately from Definition 
m Lemma 10.6] and [SJ Lemma 4.75]. 

Lemma 2.5. Let D* _ be the dual net of a polynomial lattice point set P2™, s (q 1 p). 
Then we have 

^ 2 ^ 1 wal fc (, J )=/ 1 zfkED ^ 



i=0 



otherwise. 
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3. Mean square weighted £ 2 discrepancy. In this section, we study the 
mean square weighted £2 discrepancy of scrambled polynomial lattice rules. In [5], 
Dick and Pillichshammer have derived the Walsh series expansion of the classical £2 
discrepancy. By a slight modification, we can rewrite the expression of the square 
weighted £2 discrepancy given in Proposition 1 1 . 1 1 as follows. 

Proposition 3.1. For any point set Pn_ s = {x , . . . , Xn_i} in [0, l) s and any 
sequence 7 = (7 u ) u c/ s of weights, we have 



N-l 



0#uc/ s fc„,j„eNi, u| \{o} 



where k u = (kj) j&u , l u = (Ij)jeu, r u (k u ,l u ) = l\ jeu r(k 3 J 3 ). Further, we have 
r(k,l) = r(l,k), and for non-negative integers < I < k with dyadic expansions 

k = 2" 1 - 1 H h 2 a -- 1 with 01 > ■ • ■ > a v > and I = H h 2 b -" 1 with 

b\ > • ■ • > b w > 0, we have 



KM) 



3 

1 



1 

" 2 a! + a 2 +2 
1 



3-4 a i 
1 



2 a 1 + b 1 +2 





if k = 1 = 0, 

if v = 1 and Z = 0, 

if v = 2 and I = 0, 

if v = w + 2 > 2 and 03 = 61 , . . 

ifk = l>0, 

if v = w, at ^ bi and 02 = 62, • 
otherwise. 



The next corollary provides an expression for the mean square weighted £2 dis- 
crepancy of scrambled polynomial lattice rules. 

Corollary 3.2. For a polynomial lattice point set P2m jS (q,p), we have 

E[^i, 2 » >T (iv» >i (g,p))] = E E ^ fe «>°)> 

where the expectation is taken with respect to all the possible scrambled point sets of 
p2 m ,s(QjP)- Further, we denote by (fe„,0) i/ie vector from Nq with all the components 
whose indices are not in v replaced by zero, and we have ip(k) = l/4 ai for k 6 N 
with dyadic expansion k = 2 ai ~ 1 + • • • + 2 a " _1 with a\ > ■ ■ ■ > a v > 0, ip(0) = 1 and 

Proof. Let y, y' S [0, 1) be two points obtained by applying Owen's scrambling to 
the points x,x' £ [0, 1). From Owen's lemma Lemma 13.3], we have 

E[wal fe (j/)waLV)] = 0, (3.1) 

whenever k ^ I. In the following, we denote by y n j the point obtained by applying 
Owen's scrambling to the point x n j. Using Proposition 12.31 Proposition 13.11 (|3. II) 
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and the linearity of expectation, we have 

E[£l )2miT (P 2miS ( g ,p))] 

= Y lu Y r u (k u ,l u )^ Y n E { Wal k:(yn,])™&k:j(y n >,j)] 

0^«C/, fc u ,J„ e Nj ) u| \{0} n,n'=oje« 

1 2"-l 

$^uCI s fc u eNj, u| \{0} n,n'=Oj£« 

1 1 2m_1 

Now we need to introduce the following notations. For l v = (lj)jev G N' u ', we 
define a set £>z t , as 

Si„ = {(kj)jev G N^' 1 : fo^ 1 < % < for j e t>}. 
We denote by &i v a sum of v v (ky) feu) over all fe^ £ ^z^- ^Ve have 

k„eBi v jev 

2 l i -1 

=n e r (*j.*j) 

2 ij _ 2'i-i 



n 



3 • 4^ 

1 

~ 3l«l . 2l"l+l I »li ' 

where \l v \i '■= ^2j£ V lj- Further, we introduce a so-called gain coefficient, which is 
independent of the choice of k v £ Bi v , 

1 2 m -l 

n,n'— j£v 



2H-H.I1 ^ 1; 



(fe„,o)e£>*, p 



where the last equality is appeared in the proof of Corollary 13.7]. Using these 
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notations and their values, we have 



1 1 

Z T« Z 3hA^[ Z Z r «( fc «' fc -)^r Z Il E [ walfe i^ e ^'^)] 



2 2r; 



n,n f — j£u 



V -^l V V — V 1 



31 



= E Z ^(*-.o). 

0/uCi" o 0^uCti fe^eN 1 " 

Hence, the result follows. □ 

We denote the sum in Corollary [321 by 



(3.2) 



0^uC/ a 0#-uC« fc„ e N M 

(fc„,0)GD* p 



Using the property of dual net D* _ shown in Lemma 12.51 we can derive a more 
computable form of B(q, 7). In the following, we write log 2 for the logarithm in base 
2 and we set 2^0] = . 

Lemma 3.3. Let B(q,j) be given by i3.2\) . Then we have 



5(9,7) = - E 5R + ^T E E 7«II^)' 



(3.3) 



0^«C/ S 



n=0 9^uCI s j£u 



where for x E [0, 1) we set 



4>(x) 



1 _ 2L'°S2 x i 



In particular, in case of the so-called product weights, that is, 7„ = OiPti'X?'' we 



have 



Proof. Applying Lemma l2~5l to 7), we have 

1 2 '"~ 1 

5(9,7) = E E E W ! «> )pE wal (vo)W 



2 m -l 



J_ V V ^ V 17 

2m / < / j g I u I / y 1 1 

n=0 0/uC/ 3 0^-uCuj'Gi; 



^ tp(kj)wai k:i (xn,j) 



fe,=i 
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For the innermost sum, we have by following the similar line as the proof of [TJ 
Theorem 7.3] 



^ V( fc ) wa lfe( a 



fe=i 



oo 2 ! -l 

£4/ E walfc ^ 
1 — 3 • 2 L 1 °s 2 ( ;E )J 



Then, we have 



B{q,l) = 7^ E E E II 



n=0 $^uCI a 
2 m -l 

E E 

n=0 d/nC/, 



1 — 3 • 2L 1 °S2( :E r l , 3 )J 



7^ 
3l"l 



jeu 



1 - 3 ■ 2L 1 °S2( a; ™,j)J 



which proves the first part of the lemma. 

Next in case of j v — Yijev letting 70 = 1, we have 

%7)=-En?4EEn^M 



uCI s jeu 



=-n( i +|)+^Eri[ i +7^ 

j=l n=0 j=l 

Hence, the result for the second part of the lemma follows. □ 
Since we have the following recursion in the inner sum of 

E ^Y['P( X n,j)='yirY[H x n,j)+ E ( 1 + lu ^ {r} 0(aV,j, 



n— uGI s jeu 
2 m -l s 



V>^uCI r jeu 



3 = 1 



0#uC7 r „x 



7« 



jeu 



for 1 < r < s, the computational complexity of B(q,-f) with the general weights 
is 0(2 ms ). In case of the product weights, on the other hand, the computational 
complexity of B(q,~f) reduces to 0(s2 m ). 

4. Construction of polynomial lattice rules. In this section, we first show to 
how to find a vector q by using the CBC construction algorithm. Then, we prove that 
an upper bound on B(q,-y) obtained by our algorithm converges at almost the best 
possible rate of N~ 2+s for all 5 > 0. Finally we end this section with the discussion 
about the strong tractability of our algorithm. 

We restrict each polynomial qj such that deg(gj) < m and qj ^ 0. In the following, 
we denote by R m the set of all the non-zero polynomials over F2 with degree less than 
m, i.e., 

Rm = {?£ F 2 [x] : deg(q) < m and q / 0}. 

It is clear that \R m \ = 2 m — 1. Further, we write q^ = (qi, . . . , qj) for 1 < j < s. The 
CBC construction proceeds as follows. 

ALGORITHM 4.1. For m,s € N and any sequence 7 = (p/u)uci a > we proceed as 
follows. 
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1. Choose an irreducible polynomial p G F 2 [a;] with deg(p) = m. 

2. Set ql = l. 

3. For j = 2, • ■ • , s, find q* by minimizing -B((q*_ l7 qj), 7) as a function of 
qj 6 R m where 



B((«J_ 1><& ),7)= E E E 



Vl^vCu fc„£N w 

<h v ,0)eDT , , 

v v, 1*- ,<y ), P 

The next theorem provides an upper bound on B(qj, , y) for the thus constructed 
polynomials q* for 1 < j < s, which converges at almost the best possible rate of 
f Qr a jj j > g j n ^.^g p roo f f ^jjg theorem, we use Jensen's inequality, which 
states that for any sequence (aj) of non-negative real numbers we have 



E' 



for any < A < 1 . Our proof is similar to the proof of [6j Theorem 4.4] but requires 
the slightly more complicated argument. 

Theorem 4.2. Letp £ ¥2[x] be irreducible polynomial with deg(p) = m. Suppose 
that q* £ is constructed by using Algorithm ^. 1\ Then for any j = 1, . . . , s we 
have 



B{q*,l) < 



1 



(2™ - l) 1 /^ 



-1 



2 2A -1 
2 2A -2 



l/A 



and for ^ v C Ij 



(4.1) 



E ^ A (fc„o)< 



fc„6N M 

(fe„,o) 6 £>;. 



(2 m - 1)(2 2A - 2)1*1 



(4.2) 



for 1/2 < A < 1. 

Proof. We prove the theorem by induction on j. For j = 1, we have 



fc=i 

2 m |fe 



< 



7{i} 



oc 2° — 1 



\ 



1/A 



E E ^ A ( fc ) 

»=1 k=2"- 1 
2 m \k 



/ 



1/A 



_ H 1 } j 2 a -" i ^ 1 . 2 _2a> 

1 1/A 



< 



3 
7{i} 



1 



2 2Am( 2 2A _ 2 ) 
1 



(2 m -l)(2 2A -2) 



1/A 
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for 1/2 < A < 1, where we used Jensen's inequality for the first inequality. It is clear 
that the proof of (|4.2I) is included in this calculation. Hence, the result follows for 

i = i- 

Next, let the statement of the theorem be true for some j > 2. Then it is enough 
to show that the statement is also true for the (j + l)-th component. In the following, 
we classify each subset u C according to whether u includes the component 

{j + 1}. If u includes {j + 1}, we further classify each subset v C u according to 
whether v includes the component {j + 1}. Then we have 

B((<?*, <?,-+!), 7) 

= E E 

0#tiC/ i + 1 Vl^vCu fc„GN M 

= E JE E ^-°) 

(fc„,o)en*. 

+ E 2 SSS li E E *(*-°) 

(Jt„.0)6fl*. 

+ E^^E E o) 

"C/j tiCti (fc l ,,/c J+ i)eNi l 'i+ 1 

= 5( 9 *, 7 ) + ^ + ^( gi+1 ), 



where we have defined 



and 



^=E^E E ^.o), 

(fc„,0)GL>*, 



:= E E E ^(*..*i + i,0). 



(fe„,fe J+ i)ei 



We note that 02 is a function of <7j+i, while 9i is not. 

Using Jensen's inequality and (|4.2[) . we obtain an upper bound on 9\ as 

A 



(fc„,0) G IT, 

J 

r 7uu{j +i}^ A ^ 1 
i _ i Z^ 



- V 3l«l +1 / ^ (2 m - 1)(2 2A - 2)1"! 
2»i _ i V 3l u l+ 1 V 1 _1 + V 2 2A - 2 
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In order to obtain an upper bound on 92{q* +1 ), we employ the averaging argument 
as in the proof of [HI Theorem 4.4]. Since we choose Qj+i which minimizes ^2(^+1) 
according to Algorithm 14. 1[ we have 



£(sj +1 ) < aS^TY E #fo+0 

qj+l£R m 

^ 2^ V 3M+1 ) 2^^,3+1, 



where we have defined 



fJ>v,j+i : — 2) 



r^y E E ^ A (fe„fc J+ i,o). 



(fc„,fc J+1 ,o)eD* q ,^. +i)p 
Due to the condition of (k v , fcj+i, 0) £ D? „ \ , fe« and must satisfy 
tr m (fc„) • q* + tr m (/c 3 -+i) • q j+ i = (mod p). 

If Aij+i is a multiple of 2 m , we always have tr m (fcj + i) = and the above equation 
becomes independent of qj+i- Otherwise if fcj+i is not a multiple of 2 m , we have 
tr m (fc J+ i) 7^ and the term tr m (fcj+i) - qj+i cannot be a multiple of p. Thus, we have 

00 

^•+1= E ^ A (%+i) E ^ A ( fe ») ( 4 - 3 ) 

fci+i=i fe^eN 1 " 

^^l^i+l tr^ffe^J-q^— (mod p) 

+ ^T-T E ^fe+i) E ^C^)- ( 4 - 4 ) 

2m t fe j + i tr m (feu)-«,,#0 (mod p) 

We obtained the second equality by applying the definition of the dual net and by 
classifying whether fcj+i is a multiple of p or not. In (|4.3[) , we have 

00 00 2 a — 1 

E vAfc,+i)= E E tf A te+o 

fej + l = l a=m+l fc i + 1= 2 a - 1 

2 m |fc 3 + i „ m| . 



^ ^ 2fl — m— 1 2 _ 2Aa 

=m+X 
-. 00 

9 (l-2A)a 



2 

a=m+l 
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In (|4.4j) . we have 

2^TT7 E ^ A (%+l) 



1 



kj+i—X 
m 2=-l 



2^1 E E 

o=l fc 3 + i=2«-i 



00 2 a -l 

^ E E v/fe +1 ) 



2 m - 



a=m+l k j + 1 =2" 



1 m ^ 00 

\ A 2 a ~l . 2~^Aa J- \ A t^a — l na—m—l\ 

i _ \ s 2 rn — 1 ' * 

a— 1 a— m+1 

1 m 1 00 

- w 1 - 2 ^— L_ V 

1 _ x) ^ 2" l+1 ^— 



-2Aa 



2(2™ - 1) ^ 
v y a— 1 

Hence, jJ>v,j+i is upper-bounded by 

= 1 V 2( x - 2A ) a v 



q(1-2A)o 



a— m+1 



a— m+1 



tr m {fcv ) -Qt, — (mod p) 



(?Ti OO \ 

1 Y^o a ~ 2A)a I 1 V 2 (1 ~ 2A)a l 

2(9™ _ 1) 2 m + 1 ' J I 

V_J ' a=l a=m+l J 



E ^ A ^<) 



tr m (fe t ;)"9 t ,7^0 (mod p) 



< 



1 2(!- 2A ) a 
2(2 m - 1) ^ 



1_ g 2(1 -2, )a \ ^ ^ 
a=m+l J fc^gNl"! 



^ 27^1) E 2(1 " 2A)a E 

V ' a=l fc^GNM 

In the last expression, we have 



2(i-2A)a 



0=1 



2 2A -2' 



and 



E ^) 



00 



1 M 



E^ A ( fc ) 
.fe=i 

00 2 a -l 



E E 



.o=l fc=2°- 1 
00 



.a=l 



(2 A - 2)M 
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Thus, wc obtain 

1 



< 



(2™ - 1)(2 2A - 2)M+ 1 ' 
which implies 

1 



(2 m - 1)(2 2A - 2)M+i' 
(fc„,fc J+1 )eNi"i+ 1 v A ' 

j + i 

for any v C Ij. By combining this result with (|4.2[) for ^ v C 7 S , which is supposed 
to be true by induction, we have (|4.2[) for any C i> C 
We return to obtain an upper bound on ^2(9^+1) 



uC-Ij vCu 

A . 



< 1 v P^ii±l> Vv 1 

- 2 m - 1 \ 3KI+1 Z^(2 2A -2)M +1 

uQIj vCLu 

1 V / 7uu{ J+ i} \ A 1 / 2 2A - 1 
2 m -l ^ I 3l u l +1 / 2 2A -2\2 2A -2 

Finally, by using Jensen's inequality again, we have 

= (B( q *, 1 ) + e 1 + e 2 ( q * +1 )) x 



"I 



<B x (q* n )+0 x 1+ e^ q *) 



< 



2™ - 1 ^ V3l"l/ I V 2 2A — 2 



2 2A _ lN W> 



— £ 



7uu{ J+ i}^ A f_ 1 , f2 2A -l 



tiC7j 



3M+ 1 J \ V 2 2A — 2 



7«u{ J+ i}\ A 1 /2 2A -1 



3 |«|+i ; 2 2A - 2 V2 2A - 2 



»l 



l -l ^ \?>\ u \) \ 

0#«c/, +1 \ 



2« _ 1 V3H/ I \2 2A -2 

0#uC/ j + 1 \ V 

for 1/2 < A < 1, which proves (|4.ip . Hence, the result follows. □ 

Remark 4.3. For j = 1, we can obtain the exact value for B(q*,~f) as 



6 k=i 

2 m \k 



1} 2 a_m_1 . 2 _2a 

a— m+1 



3 



Till _L 

6 2 2 " 
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We avoid this calculation in the above proof so as to prove both J^.lp and GT^D. Since 
the lower bound on the £2 discrepancy is given as hi. 2}) , this achieves the best possible 
rate of convergence. As one- dimensional polynomial lattice point set consists of the 
equi- distributed points x n = n/b m , n = 0, . . . , 2 m — 1, another QMC point sets such 
as SoboV and Niederreiter sequences constructed over F2 also give the same value. 

Here we mention about the strong tractability of our construction algorithm. Let 
us consider the inverse of the mean square weighted £2 discrepancy which is defined 
as follows 

N(s,e) = mm{N 6 N : E[£| JV)T (P JV , S )] < eE[£| ^(P , s )]}. 

We say that the mean square weighted £2 discrepancy is strongly tractable if there 
exist non-negative constants C and (3 such that 

N{s,e) < GV 3 , 

where C depends neither on e and s. 

Corollary 4.4. Assume that the weights 7 satisfy the condition 



B. 



1- A 



sup 



-1 



2*2-1 

9^ n 



1/A 



f)=£uCI s 3M 



< 00, 



for some A. Then the mean square weighted £2 discrepancy is strongly tractable. 
Proof. For empty point set Pq s-, we have 



E[£!o lT (Po >s )] = E ^1 



0^uC/ 3 j&u 

3l«l ' 



t]dtj 



= E 

0#uCI. 



For a polynomial lattice point set p2"» jS constructed by Algorithm 14. 11 we have 

1/A 



E[£ 



2,2"', 



(P*»,b)] < 



1 



(2 m - 1)VA 



E 

0/tiC7 s 



7m 
3M 



-1 



>2A 



1 



2 2A -2 



< 



( 2 m _ 1)1/A 



( 2 m _ 1)1/A 



E 

^2 



7« 

3M 



P TiA E[£| i0i7 (P 0iS )]. 



The last term is smaller than eE[£| j0i7 (P 0jS )] if N = 2 m > 1 + B*. A e~ 
result follows. □ 



Thus, the 



5. Numerical experiments. Finally, we demonstrate the performance of our 
constructed polynomial lattice rules. We focus on the case of the product weights, 
that is, 7 U = IljeiiTj' because of their importance in practice and the availability 
of the fast CBC construction algorithm using fast Fourier transform [24j[25]. Three 
choices for jj are considered here: jj = 1 (unweighted), jj = 0.9- 7 ' and jj = l/j 2 for 
j = 1, . . . ,s. 
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Table 5.1 

The mean square weighted £2 discrepancy for ■yj = I. 



m 




s — 


1 




s 


= 5 


s — 


50 


s — 


100 




SoboP 


PLR 


SoboP 


PLR 


SoboP 


PLR 


SoboP 


PLR 


4 


6.51c 


04 


6.51c 


04 


4.83c-02 


3.79c-02 


3.93c+07 


3.91c+07 


2.54c+16 


2.54c+16 


5 


1.63c 


04 


1.63e 


04 


1.45c-02 


1.37e-02 


1.96c+07 


1.94e+07 


1.27c+16 


1.27e+16 


6 


4.07e 


05 


4.07e 


05 


5.04e-03 


4.29e-03 


9.70c+06 


9.64e+06 


6.35e+15 


6.35e+15 


7 


1.02c 


05 


1.02e 


05 


1.27c-03 


1.32e-03 


4.78e+06 


4.77e+06 


3.18c+15 


3.18c+15 


8 


2.54c 


06 


2.54e 


06 


4.11C-04 


4.69e-04 


2.36e+06 


2.35e+06 


1.59e+15 


1.59e+15 


9 


6.36e 


07 


6.36e 


07 


1.21e-04 


1.38e-04 


1.17c+06 


1.16e+06 


7.94e+14 


7.94e+14 


10 


1.59c 


07 


1.59e 


07 


4.01c-05 


4.47e-05 


5.80e+05 


5.70e+05 


3.97e+14 


3.97e+14 


11 


3.97c 


08 


3.97e 


08 


1.15e-05 


1.28e-05 


2.89e+05 


2.80e+05 


1.98e+14 


1.98e+14 


12 


9.93e 


09 


9.93e 


09 


3.45e-06 


4.41e-06 


1.44e+05 


1.37e+05 


9.92c+13 


9.91e+13 


13 


2.48e 


09 


2.48e 


09 


1.17c-06 


1.39e-06 


7.17c+04 


6.69e+04 


4.96c+13 


4.95e+13 


14 


6.21c 


10 


6.21e 


10 


2.78e-07 


4.05e-07 


3.56e+04 


3.27e+04 


2.48c+13 


2.48e+13 


15 


1.55c 


10 


1.55c 


10 


7.98c-08 


1.31c-07 


1.76c+04 


1.59c+04 


1.24c+13 


1.24c+13 



Table 5.2 

The mean square weighted Ci discrepancy for 7^ = 0.9-' . 



m 






s — 


1 






s 


= 5 


s = 


50 


s — 


100 




SoboP 




PLR 


SoboP 


PLR 


SoboP 


PLR 


SoboP 


PLR 


4 


5 


86c 


04 


5 


86c 


04 


2.13c-02 


1.72c-02 


1.43c+00 


1.22c+00 


1.48c+00 


1.26c+00 


5 


1 


46e 


04 


1 


46e 


04 


6.25c-03 


5.93e-03 


6.27e-01 


5.16e-01 


6.47e-01 


5.34e-01 


6 


3 


66e 


05 


3 


66e 


05 


2.07e-03 


1.80e-03 


2.47e-01 


2.17e-01 


2.56e-01 


2.25e-01 


7 


9 


16c 


06 


9 


16c 


06 


5.25e-04 


5.41e-04 


9.81e-02 


8.85e-02 


1.02c-01 


9.19c-02 


8 


2 


29e 


00 


2 


29e 


06 


1.64e-04 


1.84e-04 


3.94e-02 


3.52e-02 


4.11C-02 


3.67e-02 


9 


5 


72e 


07 


5 


72e 


07 


4.73c-05 


5.23e-05 


1.60e-02 


1.41e-02 


1.66c-02 


1.47e-02 


10 


1 


43e 


07 


1 


43e 


07 


1.52e-05 


1.70e-05 


6.73e-03 


5.62e-03 


7.02e-03 


5.87e-03 


11 


3 


58e 


OS 


3 


58e 


08 


4.29e-06 


5.19e-06 


2.97e-03 


2.26e-03 


3.10c-03 


2.36e-03 


12 


8 


94e 


09 


8 


94e 


09 


1.25e-06 


1.58e-06 


1.25e-03 


8.90e-04 


1.31c-03 


9.33e-04 


13 


2 


24e 


09 


2 


24e 


09 


4.01c-07 


4.85e-07 


5.61c-04 


3.57e-04 


5.86e-04 


3.75e-04 


14 


5 


59e 


10 


5 


59e 


10 


9.89e-08 


1.43e-07 


2.13c-04 


1.41e-04 


2.24e-04 


1.49e-04 


15 


1 


40c 


10 


1 


40c 


10 


2.79c-08 


4.38c-08 


7.84c-05 


5.61c-05 


8.30c-05 


5.91c-05 



We compare the performance of our constructed polynomial lattice rules with 
that of the Sobol' sequences, which is one of the most well-known digital sequences 
over F2 [35j . Since the weights emphasize the relative importance of the discrepancy 
of lower dimensional projections, we use the Sobol' sequences as constructed in [12] . 
which should work as a good competitor. 

In Table l5.HI5.3i we show the values of the mean square weighted £2 discrepancy 
for the Sobol' sequences and our constructed polynomial lattice point sets, denoted 
by Sobol' and PLR respectively, with m = 4, . . . , 15 and s = 1, 5, 50, 100 and different 
choices for the weights. 

As expected from Remark 14.31 we obtain exactly the same values for both the 
rules for s = 1 and achieve the optimal rate of convergence, 2 _2m , despite of the 
choice of the weights. In case of s = 5, although the Sobol' sequences provide the 
slightly better results for large m, the values are comparable. For s = 50 and s = 100, 
we obtain almost the same values for both the rules in the unweighted case, while our 
constructed polynomial lattice rules outperform the Sobol' sequences in other cases. 
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